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Abstract 



In this paper, the equations of motion for geodesies in the neutral rotating Black Ring metric are 



derived and the separability of these equations is considered. The bulk of the paper is concerned with 
O ■ sets of solutions where the geodesic equations can be examined analytically - specifically geodesies 
confined to the axis of rotation, geodesies restricted to the equatorial plane, and geodesies that 
circle through the centre of the ring. The geodesies on the rotational axis behave like a particle in a 



^ ■ potential well, while the geodesies confined to the equatorial plane mimic those of the Schwarzschild 



metric. It is shown that it is impossible to have circular orbits that pass through the ring, but some 
numerical results are presented which suggest that it is possible to have bound orbits that circle 
through the ring. 



1. Introduction 



In recent years the study of black holes in higher dimensions has gained much greater prominence 
due primarily to the growing interest in string theory and the associated models where more than 
three spatial dimensions are considered [Hi!]. The extra degrees of freedom, due to these additional 
dimensions, have facilitated a much broader range of possibilities for gravitational phenomena. This 
is exemplified by the discovery of a black hole with a ring shaped horizon that can only exist 
when there are more than four dimensions [3]. This Black Ring solution was also the first solution 
to Einstein's equations that demonstrated black hole non-uniqueness i.e. a black hole that isn't 
described solely in terms of its mass and angular momentum. 

It has also been suggested that gravitational systems in higher dimensions may be of use in 
exploring four dimensional phenomena. There are many ideas of current research interest, such 
as brane- worlds [1] and the holographic principle [5], that rely on interpreting the effects of higher 
dimensional objects to explain apparently unrelated effects in the observable four dimensional world. 
Much of the most interesting gravitational phenomena occur within the vicinity of black holes, so 
obtaining a thorough understanding of black hole spacetimes in extended dimensions seems to be 
an essential step in developing these theories. The study of the geodesies of these black hole metrics 
gives an important insight into how matter and radiation will behave and is thus useful in developing 
higher dimensional theories. 

The general solution for a spherical stationary axisymmetric black hole of arbitrary dimension 
was derived by Myers and Perry in pj and it was shown in [7j and [8] that the geodesic equations 
of motion are separable. The separability of the equations of motion allows the radial and angular 
motion to be considered separately, giving n decoupled equations in terms of the test particle's mass, 
total energy, angular momenta in the various planes of motion, and some separation constants. Given 
that the equations can be completely decoupled, it is possible to construct an effective potential in 
the radial direction and thus determine the allowed regions of motion for different values of the 
energy and angular momenta independent of the angular variables. 

In order to express the equations of motion in a separable form it is necessary to use Boyer- 
Lindquist coordinates. These coordinates are similar to spherical polar coordinates but are ellip- 
soidal, rather than spherical, and thus the relationship to Cartesian coordinates is more complicated. 
The Kerr metric expressed in these coordinates is given by 

ds^ = -df + S + d9'^ + (r^ + a^) sin^ ^ d^^ + ^ (a sin^ ^ d0 - dt) ' (1) 

where S = + cos^ 9 and A = — 2Mr + a^. The transformations between these coordinates 
and Cartesian coordinates are given by: [H] 



t = t' 



X 



y 



yr^^+a? sin 6 cos 



VtM-o^ sin 6 sin 



— tan 



tan 



-1 



-1 



z = rcosO (2) 

where [t', r,9,(p) are the Boyer-Lindquist coordinates and a is the angular momentum per unit mass 
of the Kerr black hole. 
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In this metric it is possible to have timehke and null geodesies that orbit at constant radial 
distance with the angular coordinates varying. This gives spherical orbits with the radius deter- 
mined by the associated conserved quantities i.e. the energy, angular momenta, particle mass, and 
separation constants. An example of this sort of orbit is shown in the left hand plot of figure [H In 
|10j . Teo takes advantage of the separability of the equations of motion to calculate some examples 
of spherical photon motion for the 4D Kerr black hole. The djTiamics of the higher dimensional 
Kerr metrics are qualitatively similar, so Teo's work on spherical photon orbits could, in principle, be 
extended to higher dimensional black holes, although one would expect the extra degrees of freedom, 
afforded by the extra dimensions, to give more complicated classes of motion. 



Figure 1: The left-hand plot gives an example of the orbit traced out in three dimensions by a photon moving with 
constant radial coordinate. The right-hand plot shows the orbit when the angular momenta are chosen so that it is confined 
to a plane. 

It is also possible to choose the angular momenta so that the geodesic is confined to a single plane 
of motion at constant radius. This forces the geodesies to orbit within the plane, as in the second 
plot of figure [TJ This plot has the separation constant chosen in such a way that any geodesic, with 
^0 = 0, will orbit in the plane i.e. only will vary with time. 

The geodesies of the Kerr metric also demonstrate "pseudo-radial" motion, as shown in figure El 
This is where the test particle moves along a line of constant 9 but r and vary. In some spacetimes 
it is possible to choose the angular momentum so that (f) is also constant for all values of r, but 
the Lense-Thirring effect precludes this in the Kerr spacetime, hence the term "pseudo-radial" to 
describe these geodesies. The motion of the test particle in the left hand plot of figure [2] appears 
to describe a curve but this is because the plot is a projection of the geodesic on to a constant cj) 
cross-section. If the coordinate system is made to rotate with the test particle, so that cj) is effectively 
constant, then it becomes obvious that the test particle's motion is radial. This is shown in the right 
hand plot of figure [21 The angles 9q for which radial motion can occur are given by the equation 



where Q{6) = is the equation of motion describing the variation of 9 with time. If the energy 
E > m?, where m is the test particle mass, the geodesies are unbounded and can reach infinity, 
whilst the geodesies with E < m? are bounded from above and cannot escape to infinity. 

Given all that is known about the geodesies of the Kerr black hole it seems logical to go on and 
investigate the corresponding situation for Black Rings to see whether any of the properties of the 





(3) 
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Figure 2: An example of particle motion along a "pseudo-radial" geodesic in the Kerr metric, with constant 9 = j. The 
left hand graph plots the projection of a "pseudo-radial" geodesic on to a constant </) cross section in Cartesian coordinates. 
The right hand graph plots the same geodesic when the variation is suppressed, thus demonstrating that the geodesic is 
radial. 

Kerr geodesies are shared by those of the Black Ring. The physical properties of the neutral rotating 
Black Ring solution have been extensively studied in papers such as [H] and [3j but very little is 
known about the geodesies associated with this metric. There has been some rudimentary work in 
some papers such as [12], [13] and [H], but the geodesic calculations are very much secondary to 
the other facets under consideration. The calculations in these papers were very restricted in their 
application, so one would expect more general classes of geodesies to be more complicated. 

The Black Ring has horizon topology S"^ x S^, as opposed to the topology of the 5D Kerr 
horizon, giving greater scope for interesting classes of geodesies. This paper will explore these 
different situations as well as seeing whether any of the special classes of Kerr geodesies described 
above can be reproduced in the Black Ring metric. The toroidal nature of the Black Ring means 
that there won't be any spherical orbits, but the analogous situation, where the geodesic remains at 
a constant distance from the event horizon, will be examined. 

This paper is divided into 6 main sections. The Black Ring metric and some of its properties are 
investigated in section [2l It also describes the toroidal coordinate system for the metric and plots 
the contour lines for these coordinates. Section O derives the geodesic equations for the Black Ring 
metric and presents the conserved quantities associated with the symmetries of the metric. The 
remaining sections H] through [7] investigate some specific classes of geodesies where the equations of 
motion become separable. 

Section H] investigates the case where the geodesies are confined to the rotational axis of the 
ring. In this case the geodesies are restricted to have angular momentum in only one plane and 
their motion can broadly be divided into three separate classes. In the first case, when the angular 
momentum is zero, the geodesies pass through the origin and, in the timelike case, oscillate back 
and forth through the ring. If the angular momentum is small but non-zero, then both the null and 
timelike geodesies can be made to oscillate back and forth through the ring. In the final case, where 
the angular momentum is large, both the timelike and null geodesies are repelled by the Black Ring 
and thus can only ever be made to pass through the ring once. 

Section [5] describes geodesies that orbit at a constant radius in the equatorial plane of the ring. 
In this case the geodesies can only orbit in the same plane as that of the ring, and maintain a 
constant radius from the origin. These geodesies can also only have angular momentum in one 
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direction in order for them to remain on the plane. It is shown that there are no geodesies of this 
type in the interior equatorial plane of the ring, meaning that the only non-trivial geodesies are in 
the exterior equatorial plane. The geodesies in the exterior plane can be further sub-divided into 
null and timelike geodesies, with the only valid timelike geodesic being on the ergosurface. There 
are less constraints for the null geodesies of this type, so it is sometimes possible to have two planar 
orbits for the same angular momentum. In general the number of solutions for the null curves varies 
between zero and two, depending on the particular angular momentum of the geodesic. 

The Elth section looks at the case where the geodesies are restricted to move along circles of 
constant radius through the ring. These geodesies are constrained so that they follow contour lines 
that circle around the ring, much like a wire wrapped around a circular solenoid. In this case, the 
calculations show that there are no physically valid geodesies that execute this motion. However, 
numerical simulations of the geodesies suggest that it is possible to have bound orbits where the 
particle's radius varies as it orbits the ring. 

The final part, section [TJ calculates which classes of geodesies can perform "pseudo-radial" mo- 
tion. In the Black Ring coordinates general "pseudo-radial" lines become curves because the coor- 
dinate system is symmetric about the origin, rather than the edge of the which defines the ring. 
However, the derivation given in this section shows that this particular type of radial geodesic can 
only exist in the equatorial planes where the geodesies are straight lines. 



2. The Metric and Black Ring Coordinates 



The metric for the rotating Black Ring solution is given in terms of the toroidal coordinates 

(t, X, y, (f), i/j), where 

-oo<y<-l -l<x<l 



and the metric is given by |T5 

where 
and 



F(0 = 1 + 



Fiy) 



dy'^ ^ dx^ ^ G{x) 



G{y) G{x) F{x) 



G'(0 = (l-a(l + ^0 



C = Jxi\-iy) 



1 + A 
1 - A 



(4) 

(5) 

(6) 
(7) 



R has the dimensions of length and can be interpreted as the radius of the ring. The parameters A 
and u lie in the range 

< z/ < A < 1 (8) 

and describe the shape and rotation velocity of the ring, as described in [15]. The ring rotates in the 
ijj direction and in order for the string-like tension of the ring [T5j to be balanced by the centrifugal 
force, A has to be chosen such that 

All of the examples shown in this paper will consider Black Rings that are in equilibrium, so A will 
always satisfy the above equation. 
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In ([5]) the event horizon is given by y = — and there is a spacehke curvature singularity at 
y = — oo. The Black Ring is also rotating so it has an ergosurface where dt becomes spacelike. The 
ergosurface is given hj y = —j and, like the event horizon, has topology S"^ x S^. 

The coordinates used in ([5]) are adapted to the shape of a ring with radius R, where -0 parame- 
terises the and (x, y, (j)) parameterise the remaining space. In these toroidal coordinates fiat space 
is given by 



{x - yf 



{y' - + #T + + (1 - x'W 

y — 1 1 — 



(10) 



It isn't immediately obvious that this metric has no curvature, but a change of coordinates allows 
the metric to be transformed to something more familiar. The above metric can be transformed to 
polar coordinates using 

R^ + rf + rl ^ _ R'-rf- r| 

'(rf + ri + i?2)2 _ 4/^2^2 ^(rf + + R^y - AR^l 



This re-casts the metric in the form 

ds"^ = drl + ri#2 + drl + rldip'^ (12) 

which describes fiat space using two radial coordinates and two angular coordinates. The transfor- 
mations given in fllip give limits on the {x, y) coordinates, as ri and r2 have to remain real and 
non-negative. This is why x and y are restricted to the ranges given by (jl]). 

To understand how (x, y) span a constant ip cross-section it is useful to consider the inverse 
transformations of (JTTl) given by 

ri = R- r2 = R- (13) 

[x -y) [x- y) 

This allows the three dimensional cross-section of constant ip to be parameterised by a single radial 
coordinate and two angular coordinates. Combining the two transformations from ([T^ gives 

rfT7i = ^^^^ (14) 
x-y 

and 



tan^=^ = J^^ (15) 
ri \ I — x^ 



which are combined with (f> to span the 3D cross-section of constant ■0. 

Figure [3] uses these transformations to show the lines of constant y and x when (j) and ip are held 
constant. The ip coordinate would define the plane coming out of the page, perpendicular to the 
vertical axis. As can be seen from the plot, the lines of constant y define circles (or spheres when 
the coordinate is included) that foliate the space, with the x coordinate varying around the circle. 
The circles get bigger as y increases towards —1 with y = —1 defining the axis of rotation of the 
ring, and x = y = —1 being equivalent to r = cxd in polar coordinates. It is worth pointing out that 
the (x, y) coordinates only cover a semi-circle in one of the quadrants, depending on the values of 
and ip. Figure [3] plots the contour lines for all four quadrants to emphasise how well adapted the 
coordinates are to the shape of the Black Ring. 
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Figure 3: A two dimensional cross-section of constant (p and ijj (as well as the antipodal points </) + tt and -0 + tt) of the 
{x,y) coordinates. The red-turquoise circles (centred on the vertical axis) are lines of constant y and the blue-magenta 
circles (centred on the horizontal axis) are lines of constant x. The vertical axis on this plot corresponds to y = — 1 and 
the horizontal axis corresponds to a; = ±1, where x = +1 corresponds to the centre of the ring up to the inner edge, and 
X ^ —\ corresponds to the region from the outer edge of the ring to infinity. 
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3. Geodesic Equations and Conserved Quantities 



To obtain the geodesic equations of motion, the following lagrangian is formed 



ny) 

Fix) 



t 



CRjl + y) 
Fiv) 



^ + 



R^F{x) 
(x - yy 



l_ Gix) 
G{y) ' G(x)^F{xj 



+ 



The conjugate momenta can then be calculated using 



(16) 



(17) 



It is obvious from the form of ([5]) that there are three Killing vectors given by dt-, 9^, and 5^ 
This means that the momenta in these directions will be conserved, giving 



Pt 



-E 



F(y).^CR(l + y). 



F{x 
R^G{x) 



F{x) 



Px = 

Py 



{x - yY 

CR{l + y). C^R\l + yf 



F{x) 
R^F{x)x 
G{x){x-yy 
R^F{x)y 
G{y){x-yy 



F{x)F{y) 



G[y)R^Fix) 
^ F{y){x-yY'^ 



(18) 
(19) 
(20) 
(21) 
(22) 



where E is the energy, and I and \E' are the constants associated with the angular momenta in the 
(f) and il) directions respectively. The momenta in the x and y direction have also been included for 
future reference. 

Unfortunately, these conserved quantities aren't sufficient to allow the equations of motion to 
be separated immediately. The easiest way to check whether the equations can be separated is to 
consider the Hamilton- Jacobi equation 



a^'^Pf^Pv 



-m 



Calculating the inverse metric components from ([5]) and substituting into ( 123|1 gives 



Y{y) + X{x) 



Fix) 



'E^F{x) 



m 



(23) 



(24) 



where Y{y) and X{x) are given by 

Yiy) = 
Xix) = 



F{y) EC{l+yy 



G{y) \R F{y) 
eF{x) ^ G{x)pI 



G{y)pl 

i?2 



R^G{x) 



(25) 
(26) 



The right hand side of equation (12^ indicates that the Hamilton- Jacobi equation can't be separated 
for arbitrary values of the constants of motion but if = and m = then the function on the 
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right hand side will go to zero, allowing the equation to be separated out into terms involving only 
X and y. 

These particular values of E and m correspond to null geodesies that don't move relative to the 
space. In the case of the black ring, this means that the geodesies will rotate with the ring, but 
won't move in the x, y, or (p directions. Null geodesies with E = aren't physically realisable, but 
this does provide a way to check whether the numerical solutions to the equations of motion are 
consistent. 

Applying the variational principle to ffTBj) gives three equations of motion for the Killing directions 
as per equations ( fT8|) -(!20l). The remaining two equations of motion are calculated by varying with 
respect to x and y respectively 



H{x,y) 
H{y,x) 



J{x) 

Ay) 



e{x - yfG'{x) 

2R^G{xf 
{x-yf[ECR{l + y) + ^F{y)f 
RF{xYG{y)F{y) 



EC{1 - A) 



G'{y)F{y) 

2RG{y) ECR{1 + y) + ^>F{y) 



(27) 
(28) 



where 



C G'iOe [F{0 + F{y)]e 



G{0 2G{C? 2F(C)G(C)(C-^) 
[F{x) + F{Cm [F{Q + F{y)]r^^ 



AO 



E^X 



+ 



+ 



X 



G(C)F(x)(C-r^) 2F(C)G(r^)(C-r/) 

[EGR{1 +y) + ^Fiy)f [F{y) + F(C) - X{x - C)] e'F{0 



(29) 



2F{y) R^F{x) 



2F{x)F{y)G{y) 



G{x) 



(30) 



These two equations have been expressed in terms of the conserved quantities by substituting for 0, 
■0, and t from equations ( fT8l) -( !20l) . In certain circumstances it is also useful to use the first integral 
of motion, which is given by 



R^F{x) 



X 



+ 



e{x-yf E^F{x) {x-yf[RE{l + y)G + ^F{y)\ 



{x-yY \G{x) G{y)) R^G{x) F{y) 
where e determines the nature of the geodesies as 



F{x)Fiy)R^Giy) 



(31) 



— 1 timelike 

null 
+1 spacelike 



As previously mentioned, equation (HM isn't separable in the general case, so the following four 
sections consider some special cases where either x or y remain constant throughout the geodesic's 
motion. These specific cases give the limiting behaviour of geodesies in different parts of the space 
with varying angular momenta and energy. They can then be used to give a good idea of how the 
geodesies behave when their initial conditions are similar to any of the cases examined in sections H]- 
[71 As the initial conditions are varied away from these limiting cases, the behaviour of the geodesies 
gradually breaks down until the motion is completely dissimilar. 
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4. Geodesies Along the Rotational Axis of the Ring 



The geodesic equations, as they are presented in equations ( |T8ll28l) . are too comphcated to analyse 
straight away. In order to reduce the complexity of the problem, it is necessary to look for certain 
values of the initial conditions and conserved quantities that simplify the equations. The most 
obvious way to do this is to look for initial values of y that solve G{yQ) = 0. The reasoning behind 
this is easiest to see by multiplying (IHTl) by G{y) and then choosing y = yQ to he a. root of G{y), so 
that G{y) 0. The remaining terms are then given by 

R'n^) . . 2 ^ ix-yoy[REil + yo)C + ^Fiyo)f ^ ^ ^^2) 



(x-?/o)2 Fix)F{yo)R^ 

From this equation it is obvious that the velocity in the y direction, given by y, will be zero if the 
second term is zero. G{y) is already a fully factored cubic function with three real roots, so the 
solutions of G{y) = are given by: y^ = ±1,—-. The y coordinate is necessarily constrained such 
that —oo<y<—l, which reduces the possible values for yQ to —1 or — ^. Fortunately, yo = —1 
will cause the second term in fl32p to go to zero, provided \Ef = cEl. The line y = —1 represents the 
axis of rotation of the ring, so it is not surprising that the equations of motion become considerably 
simpler along this line. 

Ensuring that the initial value of y is zero will make sure that the geodesic doesn't move away 
from y = —1 immediately but, for the geodesic to remain on the line y = —1, y also has to be zero 
for all subsequent times. To check that y = 0, multiply fl2H|) by G{y) and substitute y = 0. This 
gives 

Fiy)G'iy)ix - yf [RE{1 + y)G + ^!F{y)f EGjx - yfjl - X)[REil + y)C + ^Fjy)] 
2R^F{x)G{y) RF{x) 
R^F{yfF{x)y (x - y) [RE{1 + y)C + ^F{y)f [2F{y) - X{x - y)] ^ ^ 



(x - yf 2R-^F{x) 

At first glance it looks like substituting y = —1 and \E' = will ensure that y = but the G{y) 
factor in the denominator of the first term causes problems because it doesn't cancel with all the 
terms in the numerator. More specifically, there will be a term of the fornJ§ 

[(l + y) + ^F(y)f 

G{y) ^ ' 

which will be indeterminate when y = —1, due to G{y) 0. To take the limit as y ^ —1, it is 
necessary to express G{y) explicitly as G{y) = {1 — y){l + y){l + uy). Expanding the numerator as 
well gives 

(l + yy 2(l + |/)vI/F(y) , ^'Fjyy 

+ -7- TT- TT- T + -7- rTZ rTZ T (35) 



(1 -y){l + y){l + uy) (1 - y){l + y){l + uy) (1 - y){l + y){l + uy) 

It is now obvious that the first two terms will go to zero in the limit as y ^ — 1 but the only way 
to ensure that the third term doesn't blow up is to define = 0. Since, y = —1 corresponds to 



^It also appears that 4" and E can be chosen to effect the same outcome for y = — ^ but, as will be seen later, 5" 
has to be set to zero if G{y) = 0. 

^The following analysis is slightly cavalier. A more detailed analysis of the singularities caused when y = — 1 is 
given in Appendix A 
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the axis of rotation of the ring, one would expect that \1/ would have to be zero because the angular 
momentum of the particle is zero when it is on the axis. This can be quickly verified by substituting 
y=-l into CT. 

Equation (!35l) also explains why y = —j; can't be used as an initial condition because, from ( |32l) . 
y can't be made to go to zero while \l/ = 0, unless y = —1. Furthermore, the line y = corresponds 
to the event horizon, so this possibility can be excluded on physical grounds. 

Having ensured that y = for all motion along the axis of rotation, it is now much simpler to 
calculate how the geodesic varies in the x direction. The x evolution is calculated by substituting 
y = and y = —1 into ( |271) and then integrating it numerically. Before doing that, it is helpful to 
calculate an effective potential for the motion along the axis to get some idea of the allowed motion. 

To calculate the effective potential substitute y = and y = —1 into (l3Ti) . Rearranging and 
expressing in terms of Px now gives 

, fFjx) _ eR'Fjx) _ E'R'FjxY ^ 

G{xY {x + lYG{x) {x + lYG{x){l-\) ^ ' 

This form of the equation can be compared with the equation of motion for a classical particle, with 
unit mass, in a one dimensional potential i.e. 

\pi + V{x) - ^ = (37) 

where £ is the total energy of the particle. In this case the effective potential V{x) can be found by 
solving for £ when = i.e. when the total energy is the same as the effective potential. Equation 
(l36ll is not quite of this form, since the equation is quadratic in E, but it is possible to construct 
an effective potential in an analogous way by setting Px = ^ and then solving for E to find two 
solutions V±{x). The motion of a particle in this potential is now possible only when E > V+{x) or 
E < V4x). 

The effective potential for (l36l) can now be calculated, giving 



Vjx) = ±. 



ejx + mi-x) _ 6(1 - A) 

\ G{x)R^F{x) F{x) ^ 



In the following it is assumed that £" > so, in this case, the only relevant potential is V{x) = V-^{x). 
The position of the turning points in the potential is given by the solution to 

^ = (39) 
ax 

The general form of this equation when R = 1 and A = Ac is 

Sx^ + Tx^ + Ux'^ + Vx + W = (40) 

where 

S = 2ez/3 (41) 

T = {Afu^ + Aeiy^-Aei^^) (42) 

U = {Afu^ + 3fiy + fu^ + 2€iy^ + 2eu-8eu'^) (43) 

V = {Aeu^ + 6fu + 2fu^ -Aeu-Afu'^) (44) 

W = {2f + 2fu^ - fu^ -3fiy + 2eu) (45) 
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In general this is a quartic equation so it is best to solve it for specific values of v and i. It reduces 
to a cubic for e = but the general solution is still too cumbersome to manipulate algebraically. 
Solving fHUj) gives the value of x for which a test particle will remain stationary. To calculate the 
minimum energy a particle can have, the solution to equation ( HOl) has to be substituted back into 
(EH]). 

J^.l. Timelike Geodesies on the Rotational Axis 




-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 -1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 



Figure 4: The left hand plot shows the effective potential for timelike geodesies with ^ = 0,1, 2, 3, 4, when u = 0.5. The 
lowest curve corresponds to £ — and £ increases with each consecutive curve. The right hand plot show how the same 
potential varies for ly = 0.1,0.3,0.5,0.7,0.9, when £ — 4. The lowest curve corresponds to = 0.9 and decreases with 
each consecutive curve, i? = 1 in both of these plots. 

Having calculated the effective potential and its turning points, it is now simple to deduce the 
shape of the geodesies on the y = —1 axis. Substituting values for i and z/ into (IHSl) gives the 
effective potentials shown in figure HI These plots show that, in general 

• the geodesies with low angular momentum penetrate further toward the origin at x = 1. 

• the potential goes from being attractive at £ = to wholly repulsive for large values of i. 

• the potential has a local minimum when u is large. 

• there is a local maximum for certain values of u. 

• when i ^ 0, the centrifugal barrier is infinite at x = 1. 

Figure [5] shows how the effective potential varies for timelike geodesies when the angular mo- 
mentum is varied between and 4 for two different values of z/. In the case of ly = ^ and R = 1, the 
potential is initially attractive for £ = but as i is increased the centrifugal barrier at a; = 1 becomes 
infinite and then widens toward x = —1. This also shifts the minimum toward x = —1 and makes 
the well shallower. This process continues until the centrifugal barrier cancels out the potential well 
completely and the potential becomes repulsive for all values of x. Although the centrifugal barrier, 
on the right of the plots in figure HI widens with £, V{—1) = 1 for all values of £ and u. This means 
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X 



Figure 5: 3D plots showing the variation of the timelike effective potential with £ and x for ly = j and = j respectively. 
In both of these plots R= 1. 



that the potential well can usually only trap particles with E < 1 but, as can be seen from the blue 
line in the right hand plot of figure IH for some values of u there is a local maximum near x = —1. 
This maximum is more apparent in the right hand plot of figure where a ridge appears at I ~ 3.5 
and X ~ —0.9. This maximum only exists for certain values of u and /, which can be determined for 
balanced rings by analysing ( l38l) when e = — 1, and A = Ac. 

The limits on £ are in general dependent upon u, so the upper limit of £ is given by 



R{T^ + V2T + 2) [l - iy{T^ + V2T + 1] 



- 32u + 3)T4 - 2iyT^ - GV^uT^ + V2{u^ - I2u + 3)T3 + {u^ - 8u + 3)T^ - l{v - If 

(46) 

where T'^ = \/2{i' — l)i> ^. The lower limit on £ is given by 

» - (47) 



Both of these equations assume that £ is positive but it is always possible for £ to be negative, in 
which case the lower and upper limits swap over and both acquire an overall minus sign. 

Having found the limits on £ it is necessary to consider the possible values of p that produce a 
local maximum in the potential. Considering the roots of ( HOl) shows that the minimum value that 
p can take, for which there can be a local maximum in the potential, is given by = |. Below 
this value it is impossible to have a local maximum near x = —1 and the potential will increase 
continuously with increasing x. Plugging u = ^ into (l46l) and (ITfl) shows that there is only one 
possible value of £ at this point given 

One other interesting feature of equation fH6|) is that the denominator becomes imaginary for 
u > 0.654. This indicates that the potential will always have a local maximum for all values of £, so 
all rings with u > 0.654 will be able to capture particles along the axis of rotation with Vmax > E > 1, 
no matter how large their angular momentum. For rings with u < 0.654, increasing £ will eventually 
smooth out the potential well and cause the potential to be continuously increasing with increasing 

X. 
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Figure 6: The plot on the left shows E — V{x) for a timelike geodesic with £ — 0, ~ 0.9 and R = 1. The right hand 
plot shows the evolution of x with r when the particles are started at a;o = —0.900 with E — 2. The initial velocity is 
chosen so that x = 0.508. 



Figure E] shows plots of E — V{x) vs x and also x vs r, where r is the affine parameter, for a 
timelike geodesic with i = 0. The left hand plot is now dependent on the energy and gives a better 
indication of how the velocity of the test particle changes as it moves along the path, particularly 
at X = +1. The allowed region of motion is now given by the area under the x axis, so the plot 
indicates that the test particle will approach x = 1, pass through the origin, and then continue out 
to infinity at x = —1. This behaviour is confirmed by the numerical simulation shown in the right 
hand plot. 

The numerical simulation has a discontinuity at x = 1 i.e. at the origin. This is a consequence of 
the coordinate system, since there is a singularity in the equations of motion at a; = 1 and y = —1. 
Naively substituting x = +1 into (1271) causes some problems because the terms which have in 
the denominator will blow up but, so long as £ = 0, it can be shown that these terms are zero when 
transformed into polar coordinates. The fact that £ has to be zero for the particle to go through 
the origin is obvious when one considers that £ measures the angular momentum around the axis 
perpendicular to ?/ = — 1, so any point on this axis (including the origin) will automatically have 
zero angular momentum. 

Unfortunately, at the origin, the transformations to polar coordinates become undefined, meaning 
that the previous analysis isn't valid and a further coordinate transformation is required. To analyse 
the behaviour of the geodesies at the origin it is necessary to transform to Cartesian coordinates, for 
which the transformations are given in Appendix A. These transformations show that for Cartesian 
coordinates, given by {zq, Zi), Zq -/^ when i; — >• at the origin. This means that the test particle is 
still moving, even though x appears to be zero, so the particle will pass through the origin and out 
into the other side of the ring. On the other side of the ring, the potential is exactly the same but 
the particle is moving in the opposite direction in the potential, so x becomes negative. 

Rearranging ( |271) to give x in terms of the other quantities and substituting x = 1, y = —1, 
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0, and y = leaved 



, + I [(1 + A) + (1 - A)]i;^ 

G(l) + 4(1 + A) 



Transforming to Cartesian coordinates and taking the appropriate limit, gives 



lim 



G{1) zo^o 



lim 
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-0 G{x) 



(49) 



The line = corresponds to the rotational axis y = —1, where zo parameterises points along this 
line and zq = zi = gives the origin. Transforming x into these coordinates and taking the limit as 
Zi —>■ gives 

hm X = -77^:7- ^ (50) 

It is obvious from this expression that x as zq but there is no requirement that io ~* 0. 
This shows that even though x — > 0, it doesn't necessarily mean that the test particle is at rest. The 
fact that X = when a; = 1 is purely an artefact of the Black Ring's toroidal coordinate system. 
Using fH9|) and fl50|) to express fHHl) in Cartesian coordinates gives 

^ = (51) 

where the second term in (1481) goes to zero because 2^0 = at the origin. This shows explicitly that 
the point at the origin is just a coordinate singularity and that the test particle passes through it 
without anything out of the ordinary happening. 

The centrifugal barrier exhibited in the potential plots of figure H] shows that it is possible to 
have low energy geodesies that oscillate back and forth along the rotational axis. When the angular 
momentum £ is non-zero, the geodesies never reach the origin at a; = 1, meaning that the particles will 
orbit in the x-(p plane, with the minimum and maximum distances away from the origin determined 
by the potential barriers on the right and left of the potential respectively. The ability of the 
geodesies to move in the (p direction means that the geodesies still pass through the centre of the 
ring but don't reach the origin. This is because the rotational axis y = —1 is actually a plane when 
the particles are allowed to move in the direction, thus allowing the particles to pass through the 
ring without going through the origin where x = 1. 

Figure [7] gives an example of this motion when low energy particles are placed within the potential 
well. The minimum of the potential is at Xmin = —0.096 and has value E = 0.642. The middle plot 
clearly exhibits periodic motion, but the period is dependent upon the amplitude. As the amplitude 
is increased the period is also increased. This is most easily seen by comparing the period of the 
green curve (with the largest amplitude) with that of the red one (with the smallest amplitude.) 
The left hand plot indicates why this happens. The unsymmetrical shape of the potential is more 
marked further away from the minimum, so only the curves with the larger amplitude will show this 
effect. 

The potential is steeper on the right hand side of the minimum potential line than it is on the 
left. This causes the maximum displacement to be greater to the left than it is to the right, meaning 



Although y is technically zero at the origin, the term in question is ^^^^ , so slightly more care has to be taken. It 
is shown in Appendix A that choosing = is equivalent to ensuring that the particle remains on the rotational 
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Figure 7: These plots show the motion of a massive particle when it is started at different distances away from the 
minimum. The initial conditions were set up so that v = 0.8, £ — 4:, and R = 1. The left hand plot shows the potential for 
the timelike geodesies. The middle plot shows the motion of massive particles when started at 0.1, 0.2, and 0.3 away from 
the minimum of the potential in red, blue, and green respectively, with the dotted black line indicating the position of the 
minimum. The right hand graph gives a specimen polar plot showing how the distance from the origin varies with (p. The 
initial conditions are the same as for the green curve in the middle plot and r ranges from to 5. In all cases x — and 
4> — initially. 



that the particle spends longer on the left hand side of the minimum potential. This gives the x 
displacement plot a slightly "bottom-heavy" appearance, with the maximum displacement being 
greater for negative x. The effect is most apparent for small displacements from the minimum, with 
the maximum displacement becoming more equal as the initial displacement is increased. 

When this is interpreted in terms of the physical motion of a test particle, it means that the 
particle orbits slowly when it is in the exterior of the ring and accelerates as it moves through the 
centre of the ring, before decelerating again on the other side. The acceleration is most marked 
when the geodesic passes close to the origin, so the particles with the highest energies will move very 
rapidly through the centre of the ring on a flat trajectory and those with lower energy will move 
through the ring on more of a curved orbit. 

The polar plot on the right hand side of figure [7] shows how the particle, corresponding to the 
green plot, moves in the x-(p plane, with the angle from the horizontal axis given by (p and the distance 
from the origin calculated using (1141) . This gives a clearer picture of the unsymmetric nature of the 
potential as each orbit is oblate with the trace precessing anti-clockwise after every revolution. If 
the trace is plotted over a longer time period then it does eventually return to its starting point. 
The polar plots corresponding to the red and blue curves (with smaller amplitudes) in the middle 
graph of figure [7] show qualitatively similar behaviour, but the precession of the orbits isn't as large, 
due to the potential becoming more asjTiimetric further away from the minimum point. 

Figure M gives a sample of the behaviour of a timelike geodesic when the angular momentum is 
large. In this case the potential is repulsive for all values of x, so a particle initially at rest will head 
off to infinity. If the initial velocity is increased then the particle can pass through the ring, with 
the minimum approach to the origin dependent on the energy of the particle. The particle will then 
go off to infinity on the other side of the ring. 

The second plot in figure [H] shows how x varies with (p i.e. the path of the test particle in the x-(p 
plane. In this case ^ is given by 

# ^ i{l + x) 

dr R^{1 - x){l + ux) ^ ' 

This shows that <p initially varies rapidly, when the particle is close to the origin, and then asymp- 
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Figure 8: The left hand plot shows the potential for £ — A, v — i, and R—\. The right hand plot shows the motion of 
the timelike particle when it is started from rest at x = 0.5 and = 0. 



totically approaches a constant as a; — — 1. This is reflected in the right hand plot of figure [HI where 
the curve levels off at (;/> ~ 2. Physically ^ — because the particle is approaching infinity and 
thus travels further and further for each interval in (p. 

4-2. Null Geodesies on the Rotational Axis 



h14 
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Figure 9: This plot shows how the effective potential for null geodesies varies for i/ = 0.1,0.3,0.5,0.7,0.9, when £ ~ A. 
The lowest curve corresponds to u = 0.9 and i' decreases with each consecutive curve, i? = 1 in both plots. 

The potential for null geodesies is shown in figure [9l with the variation of the potential for 
different values of u plotted for permissible values of x. These graphs show most of the properties 
of the null geodesies, principally: 

• geodesies of low angular momentum have a closer minimum approach to the origin. 
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• there is an infinite centrifugal barrier at x = 1. 

• only geodesies with £ = are able to pass through the origin. 

• the potential for geodesies with large i is repulsive for all x. 

• the potential can have a local maximum near x = —1 for large values of u. 




Figure 10: 3D plot showing how the effective potential for a null geodesic varies with v and x. In this plot £ = 4 and 

i?= 1. 

The variation of the potential with z/ and x is shown in figure [TU] for null geodesies. In this case 
the angular momentum (and indeed R) is purely a scale factor, as can be seen by substituting e = 
into ( l38i) . so the plot shown in figure [TO] gives the variation of the potential with u rather than £. 
This potential shares many of the same traits as that of the timelike geodesies, even exhibiting a 
local maximum near x = —1, as can be seen in figure [TO] when ly is large. The surface near x = —1 
and u ^ 1 is where the difference between the local maximum and minimum is most pronounced, 
but the width of the peak is the smallest, so it isn't very visible in figure [TUl The maximum height 
of the peak is when u = 0.653, which explains why it appears more marked at this point on the 3D 
plot. 

For the null geodesies, i and E have no independent meaning, since the test particles on null 
geodesies are massless, so it is only the ratio -| that is important. This is the reason why only one 
potential plot is given in figure When i = the potential for the null geodesies is identically zero, 
so the null particles move as if in flat space when they go directly through the centre of the ring. 

For the null geodesies, i has no bearing on whether the maximum exists (unless of course £ = 0), 
so the lower limit on u for a potential barrier to exist, is given by the solution to 

uJ - 26uJ + 36uJ - 54z/_ + 27 = (53) 
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This has four solutions but < < 1 for the equihbrium ring, so the only pertinent solution is 
given by z/_ = 0.653. Values of v less than v_ will give potentials with no potential barrier, like the 
upper curves in the right hand plot of figure [HI For values of v greater than z/„ there will always be 
a local maximum, with the position of the peak given by the solution to (!39l) . 




Figure 11: These plots show the motion of a null particle when it is started at different distances away from the minimum. 
The initial conditions were set up so that v = 0.8, £ = 4, and i? = 1. The left hand plot shows the potential for the null 
geodesies. The middle plot shows the motion of null particles when started at 0.1, 0.2, and 0.3 away from the minimum of 
the potential in red, blue, and green respectively, with the dotted black line indicating the position of the minimum. The 
right hand graph gives a specimen polar plot showing how the distance from the origin varies with (\>. The initial conditions 
are the same as for the green curve in the middle plot and r ranges from to 6. In all cases i = and = initially. 

If £ = then the null geodesies pass through the origin of the ring with the motion being almost 
identical to that of the timelike geodesic shown in figure El The major difference between the null and 
timelike geodesies, in this case, is that the null geodesies can't oscillate back and forth through the 
ring as the low energy timelike geodesies do. These null geodesies, in the toroidal coordinates, have 
the same problem at the origin as the timelike geodesies but the coordinate singularity is resolved 
in a similar manner. 

Physically, null geodesies with large E relative to i behave in a similar manner to the timelike 
geodesies with large E. When the timelike geodesies have large energy the difference between the null 
and timelike potentials is largely irrelevant, so high energy massive particles and the corresponding 
massless particles aren't affected by the black ring at large distances. As these particles approach 
the centre of the ring at a; = 1, the centrifugal barrier comes into play, meaning that the distance of 
closest approach increases with as per spherical black holes. 

If the null geodesies have small E relative to i and v > 0.653, then the curvature of the space 
allows for the null geodesies to be captured. In this case there are two orbits at a constant distance 
from the ring. These are found by solving flH^ . the largest solution giving the position of the stable 
orbit and the next largest giving the position of the unstable orbit i.e. at the local maximum of the 
potential. 

The local minimum in the potential for ly > 0.653 means that it is also possible to have null 
geodesies that oscillate through the ring. The plots of the potential and some examples of test 
particle motion are given in figure [TH The motion in this case is similar to the timelike case but the 
potential is less symmetrical, so the differences in the period of the motion are more pronounced. 
The period of the oscillations is also longer for the null geodesies than for their timelike counterparts. 

As i is increased the potential well is gradually smoothed out until the potential is always 
repulsive and takes the form given in figure [121 The plot shown in this figure shows the path of the 
null geodesic when it is moving away from the centre of the ring. The null particles in this potential 
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behave in a similar way to the massive particles shown in figure |8] but the null particles approach 
X = — 1 in a much shorter time. The null geodesies can pass through the centre of the ring but, as 
in the £ = case they will then continue off to infinity. 
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Figure 12: The left hand plot shows the potential for £ = 4, v = i, and R = \. The right hand plot shows the motion 
of the null particle when it is started at a; = 0.5 and </> = 0, with .t = initially. 



5. Planar Circular Geodesies 

To find circular orbits in the Black Ring metric it is necessary to solve the equations of motion 
so that X and y are constant for all of the motion and thus form an orbit by rotating in the '0 and cj) 
directions. In practice, this means solving flTTl) and fpSj) so that x = y = In general flTTI) and fl2H|) 
are dependent on i, y, x, and y, so finding initial values for these variables that solve y = from 
( l28l) won't guarantee that x = when they are substituted in (1271) . If x 7^ then this will cause 
X to vary, which will in turn cause y to vary. This means that both equations have to be solved 
simultaneously to find values of x and y that will give y = x = when x = y = but attempting to 
do this in general leads to intractable expressions that are of very high order. 

The simplest way to avoid this problem is to look for values of x which give x = when y = x = 
for all values of y, thus negating the need to consider both of the geodesic equations simultaneously, 
and allowing the circular orbits to be found by solving fl28|) when y = y = 0. It turns out that 
the only way to achieve this is to set x = ±1, therefore ensuring that a; = i; = 0, no matter what 
happens to y. This choice of x = ±1 confines the geodesies to the plane perpendicular to the axis of 
rotation and also forces £ to be zero, since the geodesies can't simultaneously remain on this plane 
and have angular momentum with respect to it. This equatorial plane is split into two sections with 
X = —1 being the region "outside" of the ring and x = +1 being the region "inside" the ring. 

The constraints discussed above are equivalent to demanding that the test particle be at a 
stationary point on the effective potential, with x = 0. The effective potentials for particles on the 
inner and outer equatorial planes are given in figure [13], with the derivation of the effective potentials 
given in Appendix B. The potentials plotted in figure [13] are in terms of a transformed coordinate 
z, which can be expressed in terms of y using 
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Figure 13: These plots show some sample effective potentia 
constrained to the equatorial plane. The left hand plot is the 
when 'J' = 7. The right hand plot is for the inner plane given 

i? = 1. 



Is V+(z) for red timelike and green null geodesies that are 
potential for the outer equatorial plane given by x = — 1, 
by X = +1, when = 0.1. Both plots are for v ~ ^, and 



This coordinate transformation is used to avoid the singularity caused when y ^ i-e. the 
ergosurface, due to the F{y) terms in the denominator of ( |3T1) . In terms of z, the ergosurface is at 
z = 0, with z = oo corresponding to asymptotic infinity when x = — 1 or the rotational axis when 
X = 1. The event horizon is at 2; = tanh~^ curvature singularity is then reached at 

z = — tanh~^ A. Thus the range of z is given by 

— tanh""^ A < 2; < 00 (55) 

The example effective potentials shown in figure fT3] plot V+{z) for the inner and outer equatorial 
planes, and are roughly indicative of all of the effective potentials. Varying the angular momentum 
for the X = —1 potential increases or decreases the height of the centrifugal barrier, as one would 
expect. In this case, the only circular orbits will be at the peak of the potential, and thus will be 
unstable. If the angular momentum, given by \Ef, is decreased to zero then the centrifugal barrier 
disappears and the timelike potential is strictly decreasing as z ^ — tanh~^A, whilst the null 
potential is identically zero for all values of z. This indicates that circular orbits can only exist when 
^ ^ 0. 

The sample plot for the x = +1 potential allows one to immediately deduce that the black ring 
won't support circular orbits in the inner equatorial plane. When 7^ both the null and timelike 
potentials increase with increasing z and have a centrifugal barrier at z = 00 i.e. as they approach 
the axis of rotation y = —1. When = the centrifugal barrier disappears, allowing both null and 
timelike geodesies to reach z = 00 and thus go through the origin. In this case the null potential 
is everywhere zero and the timelike potential levels off at V^{z) = |. As the angular momentum 
is increased the centrifugal barrier dominates both the timelike and null potentials, with the two 
potential curves converging rapidly as z increases. This has the effect of making the null and timelike 
potentials look identical for large \E'. 

When \E' is negative, both the potentials for the inner and outer equatorial planes are qualitatively 
similar to those given in figure [T31 The major difference is for the potential of the outer equatorial 
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plane, where the centrifugal barrier moves from being outside the ergosurface to being between the 
ergosurface and the event horizon. This is due to the frame dragging effect, whereby an incoming 
particle can make a closer approach to the event horizon if it is moving in the opposite direction to 
the rotation of the black hole. 

Having chosen the values of x and x such that the geodesies are confined to the plane, it is 
necessary to choose values of y that ensure that y = and y = i.e. the particle has to be on the 
peak of the centrifugal barrier. This ensures that the orbits will close up after each rotation. To 
find the values of y that solve y = it is sufficient to consider the first integral equation given by 
( l3Tll . Substituting ^ = 0, x = ±1 and i; = into this equation simplifies it substantially. After 
rearranging to isolate y and solving for = 0, the equation becomes 

E\l±\)G{y,) {±l-y,f[RE{l + y,)C + ^F{y,)f , „ .... 

F^) + (1 ± X)F{y,)R^ +eG(,o)-0 (56) 

Solving this equation for y^ gives values of yo that are on the effective potential line when x = ±1. 
Equation fl3U|) is a cubic in y^ due to the Giyo) coefficients, but for x = ±1 a factor of =F 1) can be 
removed. This reduces the equation to a quadratic making it much easier to analyse. Unfortunately, 
the coefficients of y are very complicated, so it is easiest to choose particular values for the conserved 
quantities and then solve the equation. 

To find the values of y for which the particle is at a turning point in the potential, it is necessary 
to simultaneously solve f^Hj) ioi y = ij = x = C = Substituting in x = ±1 leaves 

^(±1 - 2/o)'C(l - A) [RE{1 + y,)C + ^F{y,)] E^\{1 ± A) 



R{l±X)G{y,) 2 
{±l-yof[RE{l + y,)C + ^F{yQ)f \F{y,)G'{y^) , 2F{y,) - \{±l - y^) 



2R\l ± A)G(yo) 



+ 



G{yo) (±1 - yo) 



(57) 



This is a quartic equation in y but, once again, it is possible to factor out F{yo)'^, leaving a quadratic 
in y. Expanding the functions in (JSTj) and re-writing in fully factored form gives 

(1 + XyoYil - Jy){a±y^ + (3±yo + 7±) 



2i?2(i + j,2y^x -!)(! + iyy,y{l ± yo)2(±l + A) 



(5^ 



where 

a 



= 2iy{l + ufR^E^-2iy{l + uy^2{l-iy^)R'^E + iy'^'^{iy-l){u'^-Au-l) (59) 

a_ = 2u{7iy^ + u'^ + u - l)R^E'^ - 2//(l + i/)(3// - 1)^2(1 - u'^)R^E - 

z^^2(z/- l)(z/2 + 4z/- 1) (60) 

(5+ = 4z/(l + ufR^E^ - 4z/(l + uyj2{l - u^)R'^E - 2iy^^{i^ - + 3) (61) 



p_ = Aiy{iy^ + 7u'^ -u + l)R^E^ -Au{l + iyyj2{l-i^^)R^E-2iy^'^{u-l){iy^ + 3) (62) 



7+ = 2i^{l + uyR'E^~2u{l + uyj2{l-u^)R'^E~^\u-l){iy^ + 2u'-u + 2) (63) 



7_ = 2iy{-iy^ + iy^ + 9u-l)R^E^-2iy{l + iy){3-iy)^2{l-u^)R-^E + 

^2(z/-l)(iv^-2z/2-i/-2) (64) 

It is obvious from (1581) that yo = —j is always going to be a solution to this equation, so the roots 
of the quadratic part will give the non-trivial solutions to ([57]). The roots are given by 



yo = (65) 
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At this point it is worth checking that these roots can become real for a; = 1 and x = —1 because 
the imaginary results are unphysical. The condition for the roots to be real is given by 



_ 4a±-f± > (66) 

For R = 1 and x = +1 this becomes 

- 8^V(1 + iyf{u - 1)^ 2(1 + iy)E^ - 2^^2(1 - iy^)E + ^^(1 - u) >0 (67) 

The factor in front of the square brackets is always negative so, apart from the trivial solutions at 
1/ = ±1 and \1/ = 0, there can only be real solutions when the term in the square brackets is negative. 
The coefficient of the E'^ term is always positive, so the term will become negative for values of E 
between the two roots. Solving for E gives a repeated root at 



This means that the term in the square brackets will never become negative and thus there will only 
ever be a single real solution to (!58l) for y. Substituting E± into the quadratic part of (!58l) gives 

2(z/- 1)2(1 + //2/)2^2 = (69) 

which has only one solution at y = — i.e. on the event horizon. Equation ( !56ll can only be satisfied 
for null geodesies at this point, so it is a trivial solution. The solution where \1/ = is equally trivial 
because ( l58il then reduces to 

2E^iy{l + uf{l + yf = (70) 

This only has a solution at y = —1, which corresponds to asymptotic infinity. Therefore there are 
no non-trivial solutions for geodesies inside the ring at x = +1. This is in agreement with the 
conclusions drawn from the potential plot given in figure [T3l 

When X = — 1, (lUUI) is quartic in E and thus has real solutions for various values of E, and z/. 
These solutions are explored for timelike and null geodesies in the following two sub-sections. 



5.1. Circular Orbits for Timelike Geodesies 

For timelike geodesies, one of the non-trivial solutions to flFTj) has to be discounted because it 
is always less than y = —-, meaning that this solution is always inside the event horizon for all 
values of u, E, and Geodesies starting at this point can't form circular orbits because they are 
compelled to move toward the curvature singularity by virtue of the dt Killing vector being spacelike. 
This leaves two possible solutions to fl57p and fl56p respectively. 

After discounting the invalid solution curves, the circular orbits are given by the points where 
the solution curve for (l56ll intersects with that of (|57l) . This is equivalent to finding a point where 
both y and y are zero. Figure [TH shows how these solution curves vary with u for some specific 
values of E and 

The plots shown in figure [TH show that the blue curve, where y = 0, intersects the purple curve, 
where y = 0, at = 0.11385. The values of \l/ were chosen specifically so that the value of u 
where the curves intersect is the same for both of the plots. The purple curve corresponds to the 
non-trivial solution to ( 1571) . so the particular set of values used for this plot will form a circular 
orbit at y = —1.52508. The left hand plot indicates that there will only be one circular orbit when 
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Figure 14: These plots show the solution curves for the outer equatorial plane when y — in purple and red and those 
when y = in blue. In this instance the red line also indicates the ergosurface, since this is also a solution to ([57|) . The 
green line gives the position of the event horizon at — and the black line shows asymptotic infinity at y = —1. The left 
hand plot is for vj/ = 2.00 and the right hand plot is for * = -1.29. Both these plots use E = 1.20. 



\I/ > but when \1/ is negative, i.e. when the particle is moving in the opposite direction to the 
rotation of the ring, there can be two solutions, as indicated in the right hand plot. These solution 
are represented by the points where the blue curve intersects the purple and red curves. The blue 
and red curves can only intersect when the particle is rotating counter to the ring because the red 
curve represents the ergosurface, and at this point a circular orbit can only be formed if the particle 
is moving against the rotation of the ring, otherwise the frame dragging effect causes the particle to 
rotate too quickly. 

In general, the point where the blue curve intersects the purple curve will always be at the 
turning point in the blue line. This point is where both of the roots of fl56|) converge. This can be 
understood by considering the plot shown in figure [T^ In this plot the points where the curves cross 
the horizontal axis is where y = and the turning points of the curves are where y = 0. In order 
for y = y = for the same value of y, the turning point has to be where the curve intersects the 
horizontal axis i.e. the curve must have a repeated root. The plot in figure fTSl gives some example 
curves for different values of i' either side of the critical value. It verifies that there is a curve with 
a double root between 0.1125 and 0.1250, since the blue curve has it's turning point just above the 
axis and the yellow curve doesn't quite reach the axis. This agrees with the value of ly indicated by 
the intersection of the two curves shown in figure [TH 

Armed with this knowledge, the point where the curves intersect can be calculated by solely 
considering equation (|56ll . Once the trivial solution at ?/ = — 1 is factored out the remaining equation 
is a quadratic, so the double root will be where 

(Z/ - 1)2*4 ^ 2R\u - 1)[(^2 ^ g)(3^2 + + 1) + 4l/^2]^2 ^ i6i?3^z/^2(l-z/2)(S2 ^ ^ 

+ [£;4(z/4 - 24i^3 - 181/2 _ 8z/ + 1) + 2e(l + - ISu^ - bu + 1)E^ + e2(l - i^^f] = (71) 

This is a quartic in so will technically have four solutions, but only the largest and smallest roots 
are pertinent since the intermediate solutions are always for positive values of y. Solving this for \I' 
in terms of E and z/, when e = — 1 gives 



* = ±^7=^ (2Ey^+ ^{E^ -l){3u+l± 2V2sJiy{l + iy)){l + u)] (72) 
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Figure 15: This graph plots ij against y for values of v from 0.1, for the uppermost curve, to 0.15 for the lowest curve, 
in 0.0125 increments. All the curves are plotted for E — 1.2 and = 2, so that they correspond to the left hand plot of 
figure [TH It indicates how the point where the curves cross in figure [T4] is where (|56p has a double root. A similar plot is 
formed if the corresponding negative value of ^ is used. 



This equation shows how the angular momentum has to be varied for different ring geometries 
(different values of v) and different particle energies (given by E). The two solutions for ^> represent 
the circular orbits when the particle is rotating with and against the motion of the ring respectively. 
The positive \E' solution will always be outside the ergosurface but for larger p the negative \l/ solution 
can give a circular orbit within the ergosurface. 

In figure [HI the red curve represents the ergosurface, given by |/o = — ^ = — '•^"^^ \ When the 
blue line intersects this line the circular orbit is on the ergosurface. This circular orbit will always 
exist for all values of E and unlike the solution given by ( 1721) . which is complex for E <1. 

To find the value of \l/ for which the circular orbit exists, substitute U = ~\ into for e = — 1 
and a; = — 1. This allows \1/ to be expressed in terms of E and v as 

^ ^ R[i^\l~E^)+2u{l-AE^) + l-?,E^] ^^g^ 
2E^2{1 - //2) 

Figure [16] gives some examples of the permissible values of the angular momentum about the 
rotational axis of the ring, as v varies for a range of different energies. In general increases 
with the energy and also with v. The sign of \1/ remains constant as v is varied, so the trace shown 
in figure [TH] will never cross the axis. This confirms that the particle's angular momentum has to 
always be in the opposite direction to the rotation of the ring for circular orbits on the ergosurface. 

5.2. Circular Orbits for Null Geodesies 

The analysis of null geodesies is similar to the timelike case. The solutions to f[3Ul) and f[37|) 
still give four valid solution curves, two from each equation respectively, with one of the non-trivial 
solutions to (l57Il giving unphysical positive y solutions. Figure [13 shows the equivalent plots to 
figure [TH for the null geodesies. 
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Figure 16: This graph shows how 5* varies with v to produce circular orbits on the ergosurface. The value of ^ has been 
plotted for 5 different energies ranging from E = \ \.o E = The graph shows that the circular orbits of fatter rings, 
given by larger values of v, have to have higher angular momentum in order to produce a circular orbit. 




Figure 17: In both the left and right hand graphs, the red and purple lines represent the solution curves for y = and 
the blue lines represent the solution curves for y = 0. The points where circular orbits exist are given by the points of 
intersection of these lines. The green line gives the position of the event horizon and the black line shows asymptotic 
infinity at y = —1. The left hand plot has ^ = 2.92 and the right hand plot has ^ = —1.64. The constants in both plots 
have been set to R = I and E — 1.20. 
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In this figure it is immediately obvious where the curves intersect and thus the points where 
ij = y = for the null geodesies. The values of E and ^ for these plots have been chosen so that 
the circular orbits exist for the black rings with ly = 0.11385, as in figure [TH This makes it easy 
to compare the position of the circular orbits for the null geodesies with the timelike ones. The 
main difference between the two figures is encapsulated by the blue curves in the various plots. The 
red and purple curves are identical in general for null and timelike geodesies because they represent 
the second order geodesic equations. The reason that the red and purple lines are slightly different 
between figures [H] and [T7] is because the two figures use different values for which does affect the 
second order equations and thus their solution curves. 

Figure [IT] shows that the null circular orbits, that are off the ergosurface, are closer to the 
curvature singularity (at y = — oo) than the respective timelike orbits, as one might expect. The 
difference is particularly pronounced for the negative angular momentum plots, shown on the right 
of figures and [T71 where the circular orbit is outside the ergosurface for the timelike geodesies but 
inside for the null case. For any particular shape of black ring, given by fixing the value of u, the 
null circular orbits will always be closer to the curvature singularity of the ring than the respective 
timelike orbits. 

The derivation for the relationship between the various conserved quantities is similar to that 
given for the timelike geodesies. For the null circular orbits, the critical value of \1/ that allows the 
geodesies to form a circular orbit, is found by substituting e = into flTTl) and solving for To 
calculate the position of the circular orbit it is then necessary to substitute this value of into either 
( l56l) or (1571) and then solve for y. Solving ( 17T1) for null geodesies gives 



RE ( 



^ = (^20; +^[3u + l± 2^2^ p{l + I/) j (1 + z/) j (74) 

As in the analysis of the timelike case, the intermediate roots have been discarded as they are non- 
physical. The two remaining solutions represent the circular orbits when the null geodesic is orbiting 
with and against the direction of rotation of the ring respectively. The positive \l/ solution is always 
outside the ergosurface because of the frame dragging effect caused by the rotation of the ring. The 
negative solution can be either inside or outside the ergosurface, depending on the magnitude of 

The other null circular orbit is found where the blue curve in figure [17] intersects the red y = —\ 
curve. This is easily calculated by substituting y = — j; into (1561) and then solving for ^ in terms of 
E and u. Doing this gives 

$ = (75) 

thus allowing the position of the second circular orbit to be calculated as above. 

As for the timelike geodesies, there will always be a null circular orbit on the ergosurface but 
unlike the timelike circular orbits, there will always be a second solution with angular momentum 
given by ( 17^ . This means that for null geodesies on circular orbits there will always be two possible 
circular orbits for particular values of E and u. For the timelike case there will sometimes only be 
one solution for particular values of E and v, specifically when E < 1. 

There is one specific instance where the null geodesies can only form one circular orbit for given 
values of E and u. This is when the values of (1751) and ( fMl) are the same. For the null geodesies E 
and R are only scaling constants, so the angular momenta will only be degenerate for a particular 
value of ly. Equating f[75l) and fTMl) shows that a thin ring with ly = 0.04042 will have both of the 
null circular orbits on the ergosurface. 
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6. Geodesies Orbiting through the Ring 
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Figure 18: These plots show the numerical integration of the equations of motion for a timelike geodesic started at 
ix,x,y,y) = (-0.99000,3.03072,-6.01041,0), with E = 2 and £ = ^ 0. The Black Ring has radius i? = 1 and 
ly = 0.1. The left hand plot shows the orbit of the particle in the toroidal coordinates, with the right hand plot showing 
the orbit in polar coordinates. The green lines indicate the position of the event horizon. 



For geodesies that orbit through the ring at constant y, equation ( l28i) becomes 



R^F{yfF{x)y {x - yf [ECR{1 + y) + ^F{y)\ 



F{y)G'{y) , ; 



2 



+ 



eFjy) 
x-y 



G{y){x-yY 2R^F{x)G{y) 

E'F{x) [X{x -y)- 2F{y)] ^ E{x - yfCjl - A) [RCEjl + y) + ^F{y)] ^ ^ ^^^^ 



2{x-y) RF{x)G{y) 

where y and x has been ehminated using ( l3Ti) . To find possible solutions to this equation it 
is necessary to look for specific values of the constants u, e, y, E that cause all of the terms not 
involving y to go to zero. In practice this means expanding all of the terms to give a polynomial in 
since x is free to vary while y is constrained to be a constant throughout all the motion. 
In order to get a feel for the equations without having to look for general solutions it is helpful 
to look at the special case where [ECR{1 + y) + '^F{y)] = 0. This is possible in this case because y 
is being treated as a constant and all of the other terms are constants. This means that \l/ can be 
chosen so that 

ECRil + y) 

^ = — ^''^ 

Applying this to (!76|) reduces the equation to 

R^F{yfF{x)y ^ eF{yf ^ E'Fjx) [X(x - y) - 2Fiy)] 
G{y){x-yY x-y 2{x - y) 

It can be seen straight away that the only way that the terms on the right hand side can be set to 
zero is by choosing E = e = 0, which is the same constraint as was imposed in order to separate the 
equations of motion in section [3l 
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If the constants aren't constrained in any way (other than the physical constraints) then (176|) be- 
comes an eighth order polynomial in x. This unfortunately doesn't have any solutions for physically 
applicable values for z/, e, y, E. 

Figure [18] gives an example of a timelike geodesic in the exterior of the black ring, in a reference 
frame which is rotating in the -0 direction with the particle. The particle's initial angular momentum 
is carefully chosen so that it doesn't fall straight into the black hole, but it does eventually spiral 
into the ring when the integration is continued. It is possible to keep fine tuning the initial velocity, 
so that the particle stays out of the black ring longer but in the end, the particle will either spiral 
into the event horizon, or escape to asymptotic infinity at y = —1. 

The right hand plot in figure fTSl converts the orbit into polar coordinates, given in (031) and (|T5l) . 
and then plots it using 

a = rcos6 (79) 
b = rsinO (80) 

This plot gives a more intuitive picture of what is happening to the particle. As one might expect, the 
particle initially appears to be in a stable orbit but, after approximately two revolutions, the orbit 
starts to decay and then rapidly falls through the event horizon. If the initial angular momentum 
is fine-tuned further, then it is possible to have the particle orbit the ring for a significantly longer 
period with the radius varying as it orbits. Unfortunately, the orbit always seems to decay eventually. 

The orbit shown in the right hand plot of figure [18] appears to be circular but closer inspection 
shows that it is slightly elliptical. The eccentricity of the orbit increases as the energy is reduced 
until E ~ 0.8 where it is no longer possible to find a bound orbit. It would appear from the numerical 
simulations that bound orbits can be found for all values of E greater than 0.8 though. 
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Figure 19: These plots show the numerical integration of the equations of motion for a timelike geodesic in the toroidal 
coordinates. The initial conditions are set to {x,y,y) = (-0.9900000,-5.0948494,0), with £ = * = 0, for i? = 1 and 
ly = 0.1. The other initial conditions are {E,x) = (1,0.9085145). 

If the starting value of y is fine tuned further, then the orbit of the particle looks something like 
that given in figure [19] The two plots in this figure show how y varies with x and how y varies with 
X. The left hand plot gives a more detailed view of the periodic motion as the particle orbits in a 
relatively stable ellipse. Apart from the first revolution, before the particle falls into the stable orbit. 
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it would appear that the particle is moving back and forth along a line with y depending linearly on 
x, although this isn't quite true, since the trace is slightly curved. 

The right hand plot gives the phase curve for the motion of the particle. The particle slows 
down at x = ±1, with both y and x going to zero at these two points. The plot also shows how 
the motion is fairly consistent with the curves being well grouped, apart from at the beginning and 
near the end of the plotted motion, where the curve starts to diverge from the bow-tie shape. It is 
tempting to conclude from this shape that the motion in the x and y directions is given by some 
trigonometric function but this unfortunately doesn't appear to be so. The only solvable case (given 
by substituting E = and e = into equation ( pTl) ) has solutions in terms of elliptic functions, so 
it seems reasonable to assume that the solutions for the more complicated motion would also be in 
terms of these elliptic functions. 



7. Pseudo Radial Geodesies 

It was mentioned in the introduction that the Kerr metric has some "pseudo" radial geodesies, 
where the azimuthal angle remains constant throughout the motion, so that the geodesic only moves 
in the (r, 0) plane. The analogue of this for the Black Ring metric would be to find geodesies that 
move along lines of constant shown in figure [3l 

As in the Kerr metric, it is impossible to have purely radial geodesies for the Black Ring because 
there is an analogous frame dragging effect in the Black Ring metric. Combining fl20l) and f|T8|) gives 
an expression for in terms of the constants of motion 

; (x - vY'lECRil + y) - <PF(y)] 
IPF(x)G(y) 

In order for ip to be zero for all values of y 

ECR{1 + y)- ^F{y) = (82) 

Expanding this, and collecting in terms of y, gives 

ECR - ^ + {ECR - ^A)?/ = (83) 

In order for this to hold for all values of y, it would require 

ECR = ^ (84) 

ECR = (85) 

These two equations can only simultaneously be true if A = 1, in which case the event horizon 
reduces to a three-sphere, which is a rather trivial solution. This indicates that there are no radial 
geodesies for the Black Ring metric, where ip = 0. 

If ip is allowed to vary throughout the motion, then the situation becomes the opposite of that 
investigated in the "Geodesies Orbiting through the Ring" section. In this case only x is held 
constant, which means that the equation of motion fl^Tj) reduces to 



X - yf 
2R^ 



G'(x) A 



G{x) F{x) 



_ j,.^.. _ ^G{x) [F{x) + Fjy)] ^ 

^ ' 2F{x)G{y) ^ ' 



where (IHTj) has been used to eliminate the y dependence. 

Examining the solutions to equation (|86ll indicates that there are only physically consistent 
solutions for x = ±1. This class of solutions has already been examined in some detail in the Planar 
Circular Geodesies section. 
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8. Conclusions 



In this paper the equations of motion for the geodesies of the neutral rotating Black Ring metric 
were set up and numerically integrated for some special classes of solutions. The solutions can be 
broadly separated into those that are confined to the axis of rotation and those that gave circular 
orbits in the plane of the ring. It was also shown that there are no circular geodesies that orbit 
through the ring or any "pseudo-radial" geodesies. Although it was shown that there aren't any 
circular orbits through the ring, some numerical evidence was presented that bound orbits of this 
form may exist. 

The effective potential for the on-axis solutions is very similar to that for a static black hole, with 
the potential being attractive for the geodesies with zero angular momentum. In this case both the 
null and timelike particles can pass through the origin of the ring and out to infinity, or in the case 
of the timelike geodesies, oscillate back and forth. This agrees with the analogous Newtonian case 
of a massive ring when a small test particle is placed on the axis of symmetry and then displaced 
shghtly. 

Increasing the angular momentum of the geodesies causes a centrifugal barrier to appear which 
stops the geodesies from approaching the origin of the ring, as in the case of the Schwarzschild black 
hole. Even though timelike and null particles can't reach the origin it is still possible for them to 
pass through the centre of the ring. This is because the Black Ring is five dimensional, so the axis 
of rotation is actually a plane, which means the particles can go from one side of the ring to the 
other without passing through y = —1 and x — +1. The particle motion in the x-(/) plane is similar 
to a small asteroid moving in the Sun's gravitational field. The particle can either be captured and 
orbit indefinitely, or it can escape to infinity. 

The timelike potential has a local minimum, which allows for a rich array of geodesic motion 
because it is possible to have a geodesic that is in a stable orbit in the x-(f) plane near to the centre 
of the ring. The shape of the potential well is unsymmetric so the orbit is always elliptical with the 
period of the orbit depending on the initial radius: the larger the radius, the longer the period of 
oscillation. 

The effective potential for the null geodesies is very similar to that for the timelike ones when 
the angular momentum is zero, but once the angular momentum is increased the potential becomes 
totally repulsive for small u. If i/ is large enough then it is possible to create a small local minimum 
for values of z/ > 0.653. This potential is interesting, since it means that it is possible for the Black 
Ring to have light rays in stable orbits circling through it. If u is decreased, then the null geodesic 
will always go off to infinity, no matter what the angular momenta of the geodesic is. 

In the case of the planar circular orbits, the angular momentum in the (f) direction has to be 
zero in order for them to remain on the plane through the centre of the ring. This means that 
the geodesies are confined to move in only one spatial dimension. For timelike geodesies it is only 
possible to have a constant circular orbit on the ergosurface at y = —j. This requires the energy 
and angular momentum in the ip direction to be carefully chosen though. Also, this orbit only exists 
in the outer equatorial plane. It is impossible to have any circular orbits in the interior of the ring. 

For null geodesies there is always at least one solution for tp that will give a circular orbit for all 
values of u. If the angular momentum in the i/j direction is chosen to go against the rotation of the 
ring, it is also possible to have two static orbits for the same value of u. These circular orbits do 
require a certain amount of tuning because for small values of ^' it is impossible to have any circular 
orbits, no matter what the shape and size of the ring. 

The possibihty of the Black Ring metric having geodesies that orbit through the ring at constant 
y and radial geodesies of constant x was also examined but it was shown that these cannot occur, at 
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least not for these particular toroidal coordinates, where orbits of constant y describe circles. The 
numerical evidence suggests that there may be elliptical orbits through the ring for at least one value 
of but the lack of separability of the equations of motion means that it is impossible to interpret 
these orbits quantitatively. 

It would be interesting to investigate these orbits more thoroughly, to see if the motion of 
the geodesies reveals any underlying properties of the Black Ring metric that have been thus far 
overlooked. A more systematic way of doing this might be to look for regions of the space where 
the geodesies are bounded, by numerically integrating the fully specified equations of motion for 
varying initial positions. The regions of space close to the x = ±1 planes and the y = —1 axis have 
properties similar to the results presented here, so this may provide a way of estimating values for 
the conserved momenta that could give bounded geodesies for some points. 

Another related avenue of investigation has also recently opened up due to the discovery of a 
more general solution for the Black Ring [IB] , where the ring rotates in both the ip and (p directions. 
It would be interesting to examine the geodesies in this metric to see if the extremal situation where 
both of the ring's horizons are degenerate will allow the equations of motion to be separated. This 
may then give some insight into the properties of the singularly rotating Black Ring. 
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10. Appendix A: Analysis of the Singular Terms in the Equations of 
Motion 

In certain situations the equations of motion given by (fTSj) - ([20]) and (ITH) - (12^ break down, such 
as when y = —1 is substituted into (1281) . In cases such as this, certain terms become indeterminate 
in the x, y coordinates. Once these terms are isolated, they can be analysed by transforming to 
spherical polar coordinates, as given in (HM and (|T5l) . In the case of (1281) the singularities occur in 

and J— (87) 



Giy) G{y) 

In section H] it was stated that so long as \1/ = the G{ii) term would not blow up. This is more 
evident if G{y) is converted into spherical polar coordinates. Doing this gives 

where P = Vr^TW^r^cosW+R^ . 

The rotational axis, given by y = — 1 is equivalent to 9 = 0, so it is obvious that G{y)~^ oo 
as ^ — *^ because of the sin^ 9 term. Fortunately, the term will cancel out the sin^ 9 term if it is 
initially chosen to be zero. 
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Transforming the other problematic term into spherical polar coordinates gives 



y2 _ AR^r"^ - R^)r sine - {r"^ + R^ye COS Of ^^^^ 



G{y) P'^[P - p{R^ + r^)] 

Taking the limit as ^? ^ gives 

nn.^ = ^^^-^^^ (90) 

In this form it is obvious that this term is not singular, as R^ is always positive and v <1. 

There are similar problems with equation ( |27j) when the geodesies on the equatorial plane are to 
be considered. The terms in question are 

and (91) 



G{x) G{x 

Using the same process as above the terms can be transformed as follows 

±2 _ 4i?2[r^ sin ^(i?2-r2)-f cos ^(i?2 + r2)]2 
'G{x) ~ P3[F + z/(i?2_r2)] 

The equatorial plane corresponds to = |, so taking the limit gives 



(92) 



G(x) ~ (i?2 _ ^2)2(1 + ^) (93) 

This term is very similar to (lUU]) . In this case the denominator is always positive because the 
{R"^ — term is always positive. 

Transforming the other term in (l9Tll gives 

^ ^ ^ (94) 



G{x) 4i?2r2cos2^[F + z/(i?2_r2)] 

This will obviously blow up for Q = \ unless i = 0. In a similar manner to the geodesies on the 
rotational axis, i has to be zero for geodesies on the equatorial plane, since the coordinate is 
measured with respect to the x = ±1 axis. 

Equation (l92ll also causes a problem when calculating geodesies that pass through the origin. In 
the toroidal coordinates, the origin corresponds to x = 1 and y = —1, unfortunately the transfor- 
mation given in fllSp becomes undefined. This means that a different coordinate system will have to 
be used to remove the singularity at this point. A good candidate is Cartesian coordinates. 

The transformations between the Cartesian coordinates and the toroidal coordinates are given 

by 

. ±M!ZI ,96) 

y-x 
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where zq and Zi are the Cartesian coordinates on the {x, y) plane. In these coordinates 



AR'^ zo[R'^ zq + 2zoZiZi + zq^zq - z^z^ 
AR"^ zilR"^ zi - IzqZxZq - zxzx + z^z^ 
m^z^^ {yR^ -jyzo^ - vz^^ + Q) 
AR^zi^ (z/i?2 + i,zo^ + ^zi^ - Q) 



. _ -±lh 4,Q ^ ^^O^l^l -T ^0 ^0 ~ ^1 ^OJ 

X — [))() 

. _ ^^ih ^i[ih 4-1 — ^^o^i^o ^1 ^1 -r ^0 ^ij 

y — 7p 

= ^'^^ (99) 



= ^ ^ (100) 



where Q = \j[{zi — RY + -2o][(^i + + ^o]- Expressing the terms that become singular at x = 1 
and y = —1 'm Cartesian coordinates gives 



X 



2 



_ AR^[R^zq + 2zqZiZi + zq^zq - zYzpf 

G(x) ~ (z/i?2 _ j,Zo^ - UZi"^ + Q) Q3 

y2 _ AR^[R^zi - IzqZxZq - zizi + z^^zi^ 

G{y) ^ (z/i?2 + uzo^ + ^,^^2 _ g) g3 



(101) 
(102) 



In Cartesian coordinates the origin is at 2:0 = = 0, so substituting these values into the above 
equations gives 



lim 



lim 

20^0 



lim 

21^0 



lim 

21^0 





^Zq 


G{x) 


i?2(z/+ 1) 


y' ' 


Azf 


G{y)\ 


i?2(zy- 1) 



(103) 
(104) 



It is now manifest that these terms are non-singular at the origin and are dependent on zq and ii 
respectively. 



11. Appendix B: Effective Potential on the Equatorial Plane 

For geodesies confined to the equatorial planes, given by x = ±1, can be calculated from the 
first integral equation. Substituting x = and £ = in ( 1311) gives 

i?2F(x)y2 E'F{x) {x-yY[RE{l + y)C + ^F{y)f^^ ^^^^^ 



G{y){x-yY F{y) F{x)F{y)R^G{y) 

In principle, the effective potential can now be calculated but there is a problem caused by the F{y) 
terms in the denominator. These terms become singular when y = — ^ so a coordinate transformation 
is required to make sure that the effective potential is continuous across the ergosurface. The 
transformation 

z = - tanh-i (^^^t:^^ (106) 
is continuous when y ^ —j and approaches infinity as y — > — 1. 
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Expressing fllOSp in terms of z gives 



where 

K 



= KE' + LE + M (107) 



(—A — tanh z + u -\- vX tanh z) {xX + x tanh z + 1 + X tanh z) 



M 



R'^tanhz{l - A2)2(l -tanh^ z) 
{xX + xtanhz + l + Xtanhz)'^C'^{X - if 
i?2tanh;zF(x)2(l - A2)3(l + tanh z)2(A + tanh ;z) 

2(xA + xtanh2+ 1 + Atanh2)^^C(A - 1)(1 - tanh(z)) 
R^F{xy{l - A2)2(l - tanh^ ^)2(A + tanh z) 

(—A — tanh z + z/ + z/A tanh z) [xX + x tanh z + 1 + A tanh zf'^ 
(1 - A2)(l - tanh^ z) (A + tanh z)R^F{x) 
{xX + X tanh z + 1 + X tanh zY tanh z"^^ 
'i?4F(x)2(l - A2)(l -tanh^ z)2(A + tanhz) 



:i08l 



(109) 



(110) 



The effective potentials for these planar geodesies can be calculated in a similar way to those at the 
beginning of section IH so solving (11071) for E when i = gives 



V. - -"^i^ (111) 

Technically both effective potentials need to be considered, since the L term is not equal to zero, as 
was the case for the on axis geodesies considered in section 4, but in practice it is usually possible to 
consider only since V- is usually negative for all values of z. However, if \l/ < 0, is positive 
for some values of z, in which case and the portion of V_ that is positive will be considered as 
the effective potential function. 
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